The Viscosity of Meson Matter. 
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We report a calculation of the shear viscosity in a relativistic multicomponent meson gas as a 
function of temperature and chemical potentials. We approximately solve the Uehling-Uhlenbeck 
transport equation of kinetic theory, appropriate for a boson gas, with relativistic kinematics. Since 
at low temperatures the gas can be taken as mostly composed of pions, with a fraction of kaons and 
etas, we explore the region where binary elastic collisions with at least one pion are the dominant 
scattering processes. Our input meson scattering phase shifts are fits to the experimental data 
obtained from chiral perturbation theory and the Inverse Amplitude Method. Our results take the 
correct non-relativistic limit (viscosity proportional to the square root of the temperature), show a 
viscosity of order the cubed of the pion mass up to temperatures somewhat below that mass, and 
then a large increase due to kaons and etas. Our approximation may break down at even higher 
temperatures, where the viscosity follows a temperature power law with an exponent near 3. 



I. INTRODUCTION 



The conceptual framework of this paper is a hadronic 
medium at zero baryon number, and dilute enough that it 
can be considered as a meson gas, not too different from 
a perfect fluid. We have good reasons to believe this is a 
satisfactory approximation to the state of matter in the 
debris following a relativistic heavy ion collision (RHIC) 
in the laboratory, and it may also be relevant in future 
astrophysical considerations, where the relatively simple 
extension of this work to include nucleons is definitely of 
importance. The observable particle multiplicities, cor- 
relations, angular distributions, etc., measured in RHIC, 
are customarily fitted to a few parameters within hydro- 
dynamical models for the exploding gas 0. These in- 
clude typically a set of initial conditions (phase transition 
temperature after a putative quark and gluon plasma, 
QGP, reaches the confined stage, initial energy density 
and state equation, etc.). The hydrodynamical evolution 
is supposed to control the evolution of the gas through a 
chemical freeze-out temperature, after which the particle 
composition is fixed and chemical potentials become nec- 
essary, down to a thermal freeze-out temperature, where 
the hadrons suddenly abandon their equilibrium state 
within the fluid and travel unscattered to the detector. 
This sudden transition from perfect local equilibrium to 
straight particle streaming is known as the Cooper-Frye 
prescription. Whether a "geometric" (related to the fi- 
nite size of the expanding almond-shaped region) or a 
"dynamical" (related to the local expansion rate) freeze- 
out, the scattering between individual particles governs 
the mean free path which would make sensible a hydro- 
dynamical approach. 

Mesons at low transverse momentum seem to behave 
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hydrodynamically in the sense that the observed momen- 
tum distributions can be fit with hydrodynamic models 
up to a decoupling temperature 0- Still the fits are 
worse for high-p^ and observables within HBT (Hanbury, 
Brown, Twiss) interferometry and the elliptic flow coeffi- 
cient seem to start requiring viscous corrections, at least 
within the QGP 3] . A brief, pedagogical account of some 
of these issues can be found in p|. 

We expect a better approximation to be borne out 
by considering corrections to the equations of a perfect 
fluid, allowing for a smoother transition to a free particle 
stream instead of a sharp cutoff at a given temperature. 
The ideal hydrodynamical description requires the time 
scales of microscopic processes to be much faster than 
the corresponding macroscopic fluid scale. This is not 
necessarily true for a dilute pion gas at moderate tem- 
peratures. For example, in the pion mean free path is 
found to be as long as 4-5 fermi before a thermal freeze- 
out at temperatures of order 90-110 MeV (conservatively 
low), with the chemical freeze-out temperatures in the 
range 90-140 MeV. 

Thus if a gradient in the concentration of a conserved 
quantity is present in the medium, it can be smoothed 
out by pions flying randomly without colliding too often. 
This induces a need for transport coefficients in the hy- 
drodynamical equations, in particular viscosity when the 
conserved quantity is momentum. The larger the inter- 
action, the shorter the mean free path, and therefore the 
smaller the viscosity (for a dilute gas) as a consequence 
it is interesting to also document the sensitivity to the 
choice of parameterization of the pion interaction given 
in terms of diverse scattering phase shift sets. 

To approach a microscopic calculation of the transport 
coefficients one can start from rigorous quantum field the- 
ory 6] with all the generality of Green's functions coupled 
to each other through Schwinger-Dyson equations. In the 
problem at hand collisions are mostly elastic, not requir- 
ing a coupling between states with different number of 
particles, and as the mean free paths are relatively large, 
we can employ Boltzmann's molecular chaos hypothesis, 
or decorrelation between successive collisions. It is more 
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reasonable in this case to employ the statistical formula- 
tion in terms of distribution functions /, defined below 
in section [H] These also satisfy coupled equations with 
the multiparticle distribution functions (BBKGY hier- 
archy of equations, after Bogoliubov, Born, Kirkwood, 
Green and Yuan) which the Boltzmann hypothesis trun- 
cates. Quantum and relativistic effects need to be taken 
into account for a meson gas. We thus set ourselves the 
task of solving the Uehling-Uhlenbeck equation with 
relativistic kinematics. The non-relativistic formulation 
has been published in [g. The relativistic pion gas has 
been treated before by [!j at the quantum level and by 
[l3| at the classical (Boltzmann) level as the viscosity is 
concerned. We made a brief preliminary comment about 
this system in 10 . We here report the full calculation 
and include the obvious extension to a gas composed of 
pions and particles with strange quarks (kaons and etas) . 
For this purpose we include pions, kaons and etas, which 
populate the low temperature meson gas. The other pos- 
sible states, the rho, sigma and kappa mesons, appear 
through resonances in meson scattering in our calcula- 
tion. Consistent with a sort of virial expansion we 
ignore collisions between kaons and etas, and only in- 
clude those between them and the pions. This is valid 
at moderate temperatures because their density is small 
(they conform a very dilute gas up to temperatures of 
order 150 MeV or more 

II. NOTATION AND KINEMATICS. 

We gather here our conventions and notation for the 
rest of the paper, many borrowed from [12|. Since we 
will be considering meson matter at low temperatures, 
composed mainly by pions, kaons and eta mesons, latin 
indices from the beginning of the alphabet a, 6,... denote 
the particle type and will take the values a — n, K, rj in 
the formulae below. Notice that tensor quantities such as 
Tij also carry latin indices in the range 1,2,3 for Carte- 
sian coordinate labels, which will be denoted with the let- 
ters i, j, fc, .... For these tensors, a tilde Vij means their 
traceless part, as is common use in textbooks. Greek 
indices are reserved for Minkowski space and run over 
the values 0,1,2,3. The metric g a p will be taken as 
diag(H ) (77 is reserved for the viscosity). 

In the isospin limit which we employ the degeneracy 
of each species of particle in the gas is 

9* = 3 g K = 4 g v = 1 . (1) 

We will work in the approximation (related for exam- 
ple to the virial expansion in |llj|) in which the particle 
densities dN a /dV for each species satisfies 

n v » n K , n v . (2) 

Therefore we consider binary collisions between pion 
pairs or between a pion and either a kaon or an eta me- 
son. This amounts to neglecting terms of order n 2 K , nt. 



For the elastic collisions we consider, irir — > irir, irK — > 
ttK, ttt] — > nrj the particle numbers N„, Nk, N v are sep- 
arately conserved. This is a good approximation in the 
hadron gas following a relativistic collision after chemi- 
cal freeze-out. This conservation forces the introduction 
of chemical potentials collectively denoted /i a and fugac- 
ities z a = e^^ a ~ ma \ We introduce no baryon chemical 
potential as we work at zero baryon number. The corre- 
sponding elastic cross-sections for meson scattering will 
be denoted by 

d(T a b = da(ab — > ab) 

with the following momentum assignments for the initial 
and final states: 




in terms of the Mandelstam variable s = (p + pi) 2 . The 
variable t = (p — p') 2 is related to the scattering angle in 
the center of momentum frame: 

2t 

cos 8 C M = H 5- (4) 

S Pab 

which is used below in formulae l|47|l and following. The 
total momentum is P = p + pi and the total energy E. 
Since we will choose a frame where the fluid is locally 
at rest (see below) the collision needs to be taken in an 
arbitrary frame respect to which all angles and momenta 
will henceforth be referred as the collision is concerned. 
This also leads to an additional difficulty due to a double 
valuedncss in cq. I|55l) below which is not present in the 
center of mass frame usually employed in meson-meson 
scattering. The distribution functions in phase space will 
be denoted by f a , ft,, f' a , f' b , and they are shorthands for 

fa = f a {x,p;t) . 

The normalization constants for these functions (h = 1) 
are 

SO — 

9a 
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III. THEORETICAL BACKGROUND. 
A. Hydrodynamics. 

We start by writing all magnitudes for only one particle 
species. The generalization to the three types of meson 
considered is obviously additive in this subsection and 
will be understood. The energy-stress tensor for an ideal 
fluid with local velocity field U a (x) is 



T a p = -Pg a p + uU a U 



(5) 



with the enthalpy per unit of proper volume (comoving 
volume at velocity U respect to the fixed Lagrangian ref- 
erence frame) being the sum of pressure and energy den- 
sity (also per proper volume): 

oj = P + p . 

Also conserved in the fluid's evolution is the vector field 
associated to the particle density flow: 



nil a 



(6) 



where again the particle density is taken per unit of 
proper volume. The velocity field U as seen from the 
fixed reference frame can be interpreted in terms of 
the three-velocity V by introducing the 7 factor 7 = 
-1 



as 



U = >y(l,V) . 



All quantities shall by default be refered to the comoving 
or Eulerian fluid frame, where V — 0. The ideal fluid 
continuity equation is 



d t (7x7) + V(n~/V) = 



(7) 



If now we assume the fluid to be slightly out of equi- 
librium microscopically, that is, a small departure from 
the ideal fluid, the conserved quantities need to be mod- 
ified to allow for the various transport phenomena: in a 
gas particles can move from one element of the fluid to 
another at a microscopic level, carrying their charge, par- 
ticle number, energy, momentum, etc. and maximizing 
entropy, they smooth out the gradients of these quanti- 
ties in the fluid. The macroscopic description of these 
phenomena is achieved by adding transport terms to the 
conserved currents and tensors. In particular, for the 
stress-energy tensor we add a r: 



'-a/3 



-Pg a p + ojUjJfs + t, 



a/3 



(8) 



with 



T a f3 



u? = o 



since these transport phenomena act across fluid ele- 
ments, and out of the world-line of one of them. This 
implies that in the Eulerian frame, r^o = tq^ = tqo = 0. 



To first order in the velocity gradients (that is, for 
small spatial variations of the velocity field) r is 



1 



volume term 
1 



(9) 



V ij = ^(d i V j +d j V i )--d k V«d i , 



The volume viscosity £ is usually much smaller than the 
shear viscosity r), this being the reason why we concen- 
trate on the later. This has been shown for a pure pion 
gas by Davesne in 0. We expect this result to hold in 
the multicomponent extension of the theory and accept 
it henceforth. For completeness we give the expression 
for r in a general frame of reference: 

T a/ 3 = -r]s [dpU a + d a U f3 - Uplfd^Ua 

-UclPdyUp + \d^{g a& - U a U p ) 

At the hydrodynamic level of a fluid's description, r/ is 
an empirical parameter. Its value needs to be measured 
(fitted) for the different experimental conditions of the 
considered fluid. 



B. Thermodynamics. 

We collect in this subsection the relevant thermody- 
namical properties leading to the equation of state of an 
ideal multi-component Bose gas, Since we are interested 
in the leading viscosity effects, we take the thermody- 
namical quantities to be unaffected by the interactions. 
These could be corrected if interest be found in it with 
the method of the virial expansion and the physical phase 
shifts described in 0] entailing the chemical potential fi a 
for each species, the particle mass m a , the temperature 
T and inverse temperature /3, fugacity z a , and pressure 



First we give the number density for particle of species 



a: 



n a = 7T dpp 2 — 

2?r Jo e" 

, 9a 



1 



(10) 



V e -/3!>a-™a) _ 1 



The partial pressure for species a reads 



Pa = 



-g a T 

2tt 2 



dpp 2 log 

g a T 
V 



1 _ e P(M»-E) 



(11) 



log [1 



and the total pressure is simply given by the sum of the 
partial pressures 



P 



(12) 
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We will immediately drop the Bose-Einstein condensate 
term since it is only relevant at essentially zero tempera- 
ture, and keep only the integral over the state continuum 



J 



Number Density vs. Fugacity 



in both expressions. Figure ^ shows the number den- 
sity normalized to make the species-independent quantity 
n/(gm 3 ) as a function of the fugacity. 



0.2 



0.4 



0.6 



FIG. 1: Number density for a given species as a function of 
fugacity (see text) for various values of the quotient y = ^jf- . 



The internal energy density per species is then 



9a 

2tt 2 



dpp 2 



E 



(13) 



In the absence of a chemical potential it scales as T 4 . 
Corrections due to interactions (which we consider sec- 
ond order) are to be found e. g. in [ll|]. We will operate 
with the convenient shorthanded adimensional variables 



x, y, z: 



(14) 



E = m a vT 

y a = (im a = 



Z a = e P^a--m a ) 

Defining two simple functions through quadrature: 

ln(y,z) 



o VT 



t n (y,z) = 



— ( z -l e y(VT+^-l) _ I) 

x n dx 

,o (z-^Vi+^-i) - 1) 
which satisfy a recursion relation 

yz dl n {y,z) 



t r . 



(15) 
(16) 

(17) 



2n dz 

the expressions for the partial pressure and number den 
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sity become very compact: 



Pa — '-^2^2 '3/2(2/a, Z a ) 



Air 2 



-h/2{Va, Z a ) 



(18) 
(19) 



Therefore the equation of state for a one component 
gas is 



3 t 1/2 (ya,Za) 



(20) 



and summing over all species we obtain the final equa- 
tion of state for a multicomponent ideal gas in thermo- 
dynamic equilibrium 



n I \ '3/2(2/a, Z a ) 

P = 2_^ m a n a{ya,Za)- 



3*1/2(2/0, «a) 



(21) 



in which the independent variables can be taken to be T, 
z a , or T, n a by use of ifTH)!. Equation (|2~TJ) is a general- 
ization of the classical (non-relativistic) , ideal mixed gas 
state equation 



C. Kinetic Theory. 

The distribution functions satisfy an Uehling- 
Uhlenbeck equation. This is a Boltzmann-type equation 
which includes quantum Bose-Einstein statistics in the 
final state. Since we consider three particle species, we 
have a set of coupled equations: 



= c[UJ«] + c[UJ K ] + c[UJ n ] 

-TT = C[/a',/tt] 



(22) 



neglecting KK, Krj 1 r\r\ as well as inelastic interactions, 
with the collision term 



C[f a , fb] = J da ab dp lVlcl [fj[ b (l + e„/„)(l + Uib) 

-fafu(l + M'a)Q- + £bfLJ\ (23) 

Here v IO \ is the relative velocity between the colliding 
particles and in terms of the Lorentz-invariant square of 
the scattering amplitude and phase space, we have 



da. 



1 



a b 



4E(p)E lVrc 



|T| 2 dLIPS( a ;j?',pi ') (24) 



with 



dLIPS = -L?^dn{p')S(E{p)+E 1 -E , -E' 1 ) . (25) 



For pion-pion collisions, since the particles in the final 
state are undistinguishable to the strong interactions in 
the isospin limit, a factor of \ should multiply the integral 
over da^Tr- 

The collision term is annihilated by the Bose-Einstein 
distribution function which corresponds to a gas in ther- 
mal equilibrium 



foa 



c 1 



(26) 



Small departures from equilibrium are conventionally de- 
noted by 



fa — foa + Sfoa 



(27) 



The particle and energy densities can be obtained by in- 
tegrating over momenta: 

n(x) =Jdpf , p = J dpf E(p) . (28) 

The contribution to the stress-energy tensor can be 
shown to be, summing over particle species: 



(29) 



For the purpose of evaluating the shear viscosity, which 
appears at the hydrodynamic level as a tensor of first 
order in the velocity gradient, the perturbation is taken 
to be of the form 

Xa =9a% 

with V defined in J^Jl above and as a consequence of the 
contraction with a traceless tensor, only the traceless part 
of g is relevant, so we can also take 



9ij = (PiPj - ^ l]P 2 )g a {p) 



(30) 



with g(p) a scalar function, conveniently expanded in a 
polynomial base 



g a {x)=Y J B { s a) P s {x-y,z) , 



(31) 



s=0 



the polynomials P being defined in the appendix. 



IV. SOLUTION OF THE TRANSPORT 
EQUATION 

We next show how the transport equations 1|22|) can be 
simply solved near equilibrium. We describe the method 
for only one particle species and leave the generalization 
to three or more to the reader. 
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Advective Term. 



B. Collision Term. 



To simplify the left hand side in 1)221) write it as 
df df , p 



dt dt E{p) 



•v/ 



(32) 



(v = p/E(p)) with an approximately constant perturba- 
tion: 

dtf * d t fo . 

In the proper Eulerian frame we know V = but this 
does not apply to its derivatives diV which are in gen- 
eral not null. This forces us to Lorentz transform to an 
arbitrary, nearby frame in this derivation: 



fo 



r 1 



e -H,/T e {E{p)VT=V^-pV)/T _ I 



with E(p) = y/p 2 + m 2 and [i < m. Thus, in (j33), 



df _ fo 1 

dt T 1 - e -^( B (p)-^) 



P- 



dV 
dt 



(33) 



~dP 



OP 
~dt 



+ 



d\x 
df 



a - e(p) 

T 



dT 
~dt 



This can be reduced 0] by using the standard thermo- 
dynamical relations, valid near equilibrium to first order 
in the perturbation, such as the Maxwell identities, to 
the form 



2i 

dt 



fo 



T 1 



dV 

(34) 



1 dP dT 

.__ + , Wp) _„/„)_ 



and in an analogous way 



1 



fo 



1 



PiPjVij 



T l - e -/3CE(p)-*0 
p-VP + P(E(p) - ui/n)p-VT 



(35) 



Next we profit again of the vicinity to equilibrium and 
employ the equation of continuity J7J) and state equa- 
tion l(2T!|) , while we ignore all terms proportional to VT, 
VP, V • V since they influence the calculation of the vol- 
ume viscosity or thermal conductivity, but not the shear 
viscosity upon comparison with 10 to obtain the final 
expression for the advective term: 



fo 



# = 

dt TE 1 



1 



PiPj 



v 13 



(36) 



The integrand of any of the collision functionals 
C[f a , fb] contains the product of Bose-Einstein distribu- 
tions with enhanced final state phase space: 

F = fU[ b ^Hafa){lHbhb)-fah b {l+U' a ){lH b f[ b ) ■ 

(37) 

For perturbations near local equilibrium, the distribution 
function can be written as in eq. (|27|l and this entails for 
the integrand 

F = F Q + SF . 

The collision functional evaluated on the Bose-Einstein 
equilibrium distribution is naturally zero: 



C[f0aj0b] =0 

After some simple algebra we can show 



(38) 



(39) 



5F _ f Oaf Qlbf Oaf 316 c 0(E-n. a -n h ) /^t^M _ g-/3(^-/i)\j 



A[A] = X' a +X[ b -X a -X 



lb 



and it finally follows 



C[fa,fb] = ( da ab dp 1 v rc ie l3E 

Tz a z b J 

f0af01bf0af01bA[ X (l - ^ {E - (ji))] 



(40) 



C. Expression for the Viscosity. 

Upon substitution of (|31|l on 14011 . and neglecting 5f 
against fo on the advective term 1)36(1 the linearized trans- 
port equation takes the form 



A^+AK+A^ A 12 A 13 
A 21 A 22 
-43i A 33 



BE 




(41) 



Once the integrals in A's (this is the largest computa- 
tion) and C's have been performed, the matrix system 
can be solved and the B coefficients in 1)31)1 provide the 
first approximation to the function g(p). The series con- 
verges very fast as shown in Q and therefore we keep 
only this order in the present work. This can then be 
substituted in the viscosity expressions below which can 
then be easily integrated on a computer. 



10T 



dp' , , , 

~~j^rPiPjfoa9a 



(42) 
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(compare with expression 4.19 in reference [14)1, The 
angular integrals can be trivially performed, reducing the 
expression to 

I 



which in terms of dimensionless variables (those of g can 
be read off eq. (l2*7|) ) gives 



Employing (|31|l we can write 

27r °° ® f 00 

V= -T^yZY.B^-^ dxW 5/2 (x a ;y a ,z a )P s {z a ;x a ,y a ) (45) 

i5J a s=o * a J ° 



V = -^t B ^f A °^ (46) 



and using the orthogonality relation (IB2 



V. MESON SCATTERING. 

With our kinematical choice of variables (see section 
ITT)) we can express the square of the scattering amplitude 
for 7T-7T collisions as 



the approximation studied in [l5j corresponds to 
taking a simple resonance saturation parametriza- 
tion for the isoscalar and isovector phase shifts: 



\T\ 



1 



|Tol 



sin 2 <5 00 (s) + 27 1 + 



-3|T!| 2 
2t 



5|T 2 | 2 ) 



32tt 

Ptttt 



(47) 



spin 



sin 2 Six (s) + 5 sin 2 62 



s oo(E) = - + arctan I - 

- , 7r ( E — m, 

°u{ E ) = - + arctan I p ^ 



with momentum dependent cr, width r o 
and mass m a = 5.8 m n , and p width 



(50) 



= 2.06 p 



where the isospin 0,1 or 2 has been averaged and only the 
lowest partial waves are kept, since they provide most of 
the interaction at low energies. For pion/kaon collisions 
we can also average over isospin to obtain: 



r p (p) = 0.095p 



p/mn 



1 + (p/m p ) ; 



IT 



1 



KK I 



(2|T 3 



1/2 1 



16tt 



2 sin 2 <53 + sin 2 <5i + 9 sin 2 5 1 1 ( 1 



and finally for 77- 7r scattering 



4|T 3/2 | 2 ) 
2i \ 



=(48) 



and mass m p = 5.53 whereas the scalar isoten- 
sor phase shift is simply parameterized as a straight 
line 



/167rV 
\P^vJ 



sin 2 S 



10 



(49) 



ho = -0.12p/m 7T . 



(51) 



2. 



To connect with the non-relativistic calculation, we 
will also shortly employ a (totally unrealistic at 
moderate temperatures above say a few MeV) pion 
scattering amplitude based on the low energy the- 
orem of Weinberg. This is 



Notice the partial wave normalization gives (16) for dis- 
tinguishable and (32) for identical particles in agreement 
with 14]. If only a pion gas is considered, the relevant 
phase shifts are the J = 0, I = 0,2 and I = 1, J = 1 
which dominate the cross section at low energy. We next 
study the pion, and then the meson gas, in various ap- 
proximations which are 



23 m£ 

3 fn 



(52) 



3. 



Next, through our calculation we profit from the 
IAM (the Inverse Amplitude Method 16]) fit to the 
pion scattering phase shifts. Whether in the pion 
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gas or in the mixed pion, kaon and eta gas, chiral 
perturbation theory provides an expansion of the 
scattering amplitudes at low energy: 

tu = tf)+tf} + ... (53) 

where the order of a term in the expansion counts 
the powers of m„ or momentum, or equivalcntly 
inverse powers of f n . The single channel inverse 
amplitude method constructs a model amplitude 
based on this low energy expansion which, incor- 
porating exact (as opposed to perturbative order 
by order) unitarity above the two particle threshold 
and below inelastic thresholds, allows to extend chi- 
ral perturbation theory providing satisfactory fits 
to the scattering amplitudes up to around 1 GeV. 
The one-channel IAM amplitude is 

,(2) 

*« = ( 54 ) 

i ~ 7127 
ij 

and can be understood as the [1,1] Pade approxi- 
mant corresponding to the Taylor series l|53l) . The 
first order follows from Weinberg's low energy the- 
orem. The second order includes chiral perturba- 
tion theory meson loops and order p 4 counterterms. 
These have a series of coefficients, usually denoted 
li after the work of Gasser and Leutwyler [13 which 
can at the moment only approximately be com- 
puted theoretically 01- These coefficients are in 
practice fitted to the pion scattering amplitudes or 
other observables. Since this rational approxima- 
tion to the scattering amplitudes presents in gen- 
eral poles (in particular the a and p mesons are 
clearly visible) it has been applied recently to the 
behaviour of resonances in a thermal pion bath |l9j| . 
Here we employ the fit of that paper to the scat- 
tering amplitudes at zero temperature, that is, we 
ignore the possible effect of the thermal bath on the 
parameters entering the fits (^'s, /„• and to^ which 
are adjusted to their physical value) as a higher or- 
der effect. The Us employed in our SU{2) fits to 
the pion scattering amplitudes are given in table [I] 

TABLE I: Values of the SU(2) chiral perturbation theory 
parameters employed in the IAM fit to the pion scattering 
amplitudes (input to this calculation of the viscosity). 



7i 


-0.27 


h 


5.56 


h 


3.4 


h 


4.3 



4. Finally, the full calculation including kaons and 
etas requires a parameterization of the elastic 
SU (3) phase shifts. This is again provided by chiral 
perturbation theory now with the strange quark in- 
corporated [2(j and extended to higher energies via 
a unitarization method, now the coupled-channel 
IAM. We use the phase shifts obtained in [2l| with 
the L chiral perturbation theory parameters shown 
in table |n| 

TABLE II: Values of the 5(7(3) chiral perturbation theory 
parameters employed in the IAM fit to the meson scattering 
amplitudes (input to this calculation of the viscosity). 





0.59 


L 5 


1.8 


L 2 


1.18 


L 6 


0.006 


L, 


-2.93 


L 7 


-0.12 


U 


0.2 


L 8 


0.78 



The use of the IAM is of course not necessary to 
obtain a good description of the transport coeffi- 
cients, but it is a very convenient parameterization 
of the experimental data, providing outstanding fits 
to the complete set of meson scattering channels 
described in [2l| . 

We finally quote the meson masses employed in this pa- 
per: = 139.57 MeV, m K = 493.677 MeV, m v = 
547.30 MeV. 



VI. NUMERICAL EVALUATION 

A number of integrals need to be numerically evalu- 
ated to complete the calculation. Those of the type l|B2|) 
needed to compute the coefficients C in equation (|41ll as 
well as the viscosity in (|46|l are simple one-dimensional 
integrals. The Bose-Einstein factors ensure very good 
convergence at high momentum. They are evaluated in a 
one-dimensional grid with 2000 points by an open trape- 
zoidal rule in essentially no computer time. More compli- 
cated are the integrals stemming from the collision term 
of the Uehling-Uhlenbeck equation needed for the A a b 
coefficients in <|41ll . Since we have already chosen the 
system of reference where all magnitudes are expressed 
in the comoving fluid frame, we have no freedom to study 
the collision in the center of mass, so we employ arbitrary 
momenta in these integrals. The nominal variables are 
p, pi, and fl(p'), while p\ ' is fixed by momentum con- 
servation and \p'\ by energy conservation. The later has 
a somewhat messy expression, in terms of the angle (3 
between p ' and P: 
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/, P 2 -m 2 +m 2 \ n a , /Z P 2 -m 2 +m 2 \ 2 7~ "T P 2 cos 2 ,3\ / 2 (B 2 -P 2 +mJ-m^ \ 

(! 1* «-JPcos/3±y(l P 2 cos 2 /? -4(1- ^ a j (m? 5^ 

cos/3 = — (pi s'm9(pi) sin#(p ') cos0(p") + cos6(p')(p + pi cosd(pi))) (55) 



Two solutions of this quadratic energy equation are possi- 
ble and need to be summed over in the integrand for cer- 
tain kinematical configurations. Of the eight remaining 
integrals, three are trivial. Invariance of the collision un- 
der three-dimensional rotations allow us to perform the 
angular integrals associated to, say, p and refer all angles 
to the p-axis. Then there is still an axial symmetry of the 
other three particles around this axis which allows to per- 
form the azimuthal integral over <j>(pi). The remaining 
five variables are then \pi\, 6 (pi), 0(p '). The re- 
sulting integral is performed numerically with the Vegas 
[2^| random-point algorithm. Again the Bose-Einstein 
factors concentrate the integrand in a compact set and 
convergence is very fast. For precision around 1/1000 
and x 2 ~ 1 it is enough to start with 2000 points and 
double that number 5 or 6 times, while evaluating the 
integral 10 times for each fixed number of points. 

VII. PION GAS VISCOSITY 

We first evaluate the viscosity with a constant ampli- 
tude The result is plotted in figure |2 Since the 
cross section is independent of energy, an increase in tem- 
perature simply populates states with faster pions which 
transfer momentum more efficiently. Thus, the viscosity 
grows out of control in an unrealistic manner. But this 
calculation is useful to check the non-relativistic limit re- 
ported in [8j, which used precisely this interaction. The 
behaviour at low temperature (and this is common to all 
our calculations) is governed by a non- analytic 77 cx VT 
behaviour. To see it, simply remember that for a hard- 
sphere classical gas, the mean free path is inversely pro- 
portional to the cross section and density; 

\/2na 

and calculating the momentum transfered by random 
flight of the gas molecules one can obtain 

1 

i] = —nmv v a 
3 

in terms of the rms. velocity. This shows that the 
viscosity is inversely proportional to the cross section 
and upon employing the equal partitioning of energy 
mv 2 /2 ~ 3fcT/2 we see that the viscosity is also pro- 
portional to the root of the temperature, which provides 
a convenient check. 



Shear Viscosity for Hard Sphere Pion Gas. 




FIG. 2: Shear viscosity of the pion gas with a constant scat- 
tering amplitude (from Weinberg's theorem) . Since the inter- 
action does not grow with the pion momentum, the viscosity 
is unacceptably large even for somewhat low temperatures. 
But this is used to check the low temperature limit. 



Next we turn to some more realistic pion interactions. 
These are provided by the simple analytical fit to the pion 
phase shifts from H50|l and by the SU(2) or SU (3) inverse 
amplitude method. The results are quite consistent and 
plotted in figures El 0] and \5\ respectively. The difference 
between them gives us an idea of the sensitivity of the vis- 
cosity to the employed phase shifts, since all sets of phase 
shifts are reasonable. To what precision these scattering 
phase shifts are known is an ongoing debate |23| and if 
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Shear Viscosity, Pion Gas (phase shifts from Welke et al.) 
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T(GeV) 



FIG. 3: Shear viscosity of the pion gas from the simple ana- 
lytical phase shifts (15011 from Welke et al. 

a future determination pinned them down with greater 
accuracy a much better prediction for the viscosity could 
be made, since the parameterization used for the phase 
shifts seems to be one of the largest uncertainty sources 
in the present (already realistic) computation. 



VIII. FULL PION, KAON, ETA GAS 
VISCOSITY. 

Finally we turn to a gas including kaon and eta mesons. 
We first introduce either type of particle separately and 
plot it in fig. |SJ The effect of kaons is much more dra- 
matic than the effect of eta mesons. By comparing with 
fig. which was calculated with the same pion phase 
shifts we can see that already at 100 MeV the addition 
of the kaons alone give a viscosity way bigger than present 
in the pion gas. The sensitivity of this calculation on the 
pion fugacity (density) is also very large. The combined 
effect of adding both kaons and etas to the gas is finally 
plotted in fig. [7| Of course, in a relativistic heavy ion col- 
lision we expect the chemical freeze-out of heavier mesons 



Shear Viscosity, Pion Gas (phase shifts from SU(2) 1AM). 

0.018 - 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 
T(GeV) 



FIG. 4: Shear viscosity of the pion gas from the Inverse 
Amplitude Method and 5*17(2) chiral perturbation theory. 



to occur before and therefore we need to introduce chem- 
ical potentials for all species, those corresponding to the 
heavier mesons being larger than for lighter mesons. 

In natural units the viscosity has dimensions of an en- 
ergy cubed, and indeed at moderate to high tempera- 
tures, the viscosity approximately follows a power law 
with an exponent near (and slightly above) 3, which is 
suggesting the viscosity is dominated by the highest en- 
ergy pions where the mass and chemical potential scales 
are less important than the momentum scale. The low 
temperature behaviour of the viscosity is plotted in fig- 
ure[S]where we observe how at low temperature the effect 
of adding the more massive particles is to decrease the 
plateau in which the viscosity is approximately indepen- 
dent of the temperature. 

The resulting viscosities for various temperatures and 
chemical potentials are tabulated in table ITTT1 and convey 
our final results. 
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Shear Viscosity for Pion Gas with SU(3) Phase Shifts 



Shear Viscosity for Pi/K (upper) or Pi/Eta gas (lower panel] 
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FIG. 5: Shear viscosity of the pion gas from the SU (3) Inverse 
Amplitude Method phase shifts. 



FIG. 6: Separate effect of adding K or r\ mesons to the pion 
gas, with their chemical potential vanishing. The kaons are 
by far more important (partly due to their multiplicity). 



IX. COMPARISON WITH OTHER 
APPROACHES AND DISCUSSION 

The viscosity of a pion gas has already been treated 
by D. Davesne in |9| as a quantum relativistic system. 
In this pioneering work to solve the relativistic Uchling- 
Uhlenbeck equation the pion interaction was modelled 
following ref. [l5|. which corresponds to the phase shifts 
(|5U|I above. If we consider only a pion gas and employ 
the same set of phase shifts then we can approximately 
reproduce this results. Furthermore, we document the 
sensitivity of the viscosity to the parameterization of the 
phase shifts, which was not treated in this reference. We 
have also streamlined the numerical solution of the trans- 
port equation, by employing a new family of orthogonal 
polynomials which allows to systematically extract bet- 
ter approximations if so wished, and performing a Mon- 
tecarlo evaluation of the collision multidimensional in- 
tegral, whereas the more analytical treatment in [|J is 
also somewhat more obscure. It is also interesting that 
a relaxation time estimation of the Boltzmann equation 
(with no quantum corrections) permits to approximate 



the shear viscosity of a pion/kaon gas in |14l l24j. The 
order of magnitude and qualitative behaviour as a func- 
tion of temperature are correct. We obtain somewhat 
larger results which are not unexpected as quantum cor- 
rections in a Bose gas may tend to decrease the cross 
section for scattering to initially unpopulated states, in- 
creasing the viscosity, and furthermore some difference is 
expected due to our use of the SU (3) phase shifts. 

To summarize, we have presented a systematic calcula- 
tion of the shear viscosity in meson matter at moderate 
temperatures. We have found the viscosity to behave 
as expected from a non-relativistic gas point of view at 
very low temperatures, to stabilize and even decrease (de- 
pending on the chemical potential) at small temperatures 
because of the larger cross section at increasing energies 
(the decrease would be a typical effect of a Goldstone bo- 
son gas) and to follow a positive power law at moderate 
to high temperatures. At high T (above 150 MeV) our 
approach should be less reliable because we are employ- 
ing scattering phase shifts parameterized up to momenta 
of 1 GeV and with sizeable temperatures, states with 
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FIG. 7: The shear viscosity for the full meson (n, K, if) gas 
in linear (upper) and semilogarithmic (lower) panels. At very 
low 2 and T the reader can appreciate how the non-relativistic 
behaviour (r/ oc VT) is recovered. The chemical potential for 
both K and rj has been set to zero. 



higher momentum start being populated. Eventually one 
reaches the phase transition temperature, and any results 
obtained from within the chirally broken phase (as built- 
in in our use of meson fields and chiral perturbation the- 
ory) is simply not appropriate. Extensions of this work 
to include nucleons or to evaluate other interesting trans- 
port coefficients are now straight-forward. 
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TABLE III: Shear viscosity at various values of chemical 
potentials and temperature. Units are MeV (T and /j, a ) and 
(100 MeV) 3 (rj s ). The error due to the Montecarlo evaluation 
of the 5-D integral is (5) on the last significant digit of the 
viscosity. This calculation employs the Inverse Amplitude 
Method phase shifts within SU(3). 
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APPENDIX A: LINEARIZED TRANSPORT 
EQUATION COEFFICIENTS 

In this appendix we provide the matrix elements neces- 
sary for the first order Chapman-Enskog solution of the 
transport equation. The right-hand side of the system 
gU is 




/ ^5/2(2/71-1 z tt) \ 

Anm%, 1 I \ 
-^ L h/2{VK,Z K ) 



(Al) 



The left-hand side matrix elements can be given by the 
formula 



(A2) 



Z a Zb 



daabV^idpdpxe 



/oa/ofcl/oa/ofcl 



(i-M^)(W-iWp«W 

in terms of a tensor (A. a b)kl which takes the values 

(A w u ) w =p' k p\l-e-^ E '-^) -p k p l {l-e-^ E ^-^ ) ) 

+p'iVi'(l - e-^-^)-p k lP [(l - e -^-^) 
(A K ii) H = (A f '„) fei = 

(A ia ) w = piVi'Cl - e-^-^)) -plp\{l - e-^-^) 
(Ai 3 )« =ft.Vi(l " e"^-^) - e^ 1 "^) 

(A 22 ) fc; =p'V ( (l ~ e^^'^'j-pVtl - e -0(^Cp)-Mr)) 

(A 2 i)fei = (Asijfej = 
_ e -/3(^-^)) -^(1 - e -0(^-^)) 

(A 33 )fci = pV(l - e-^ E '-^) - p k p l (l - e^™-^') 



APPENDIX B: ORTHOGONAL POLYNOMIALS. 

In solving the Uehling-Uhlenbeck equation with rela- 
tivistic kinematics and quantum statistics, one needs to 
integrate over the measure 



d/j, r (x; y, z) — w r (x; y, z)dx = 
x r dx 

A /T _ p^( z -i e y(v^-i) _ 1) 



(Bl) 



With the variable x and parameters y, z defined in 114|) 
above, with ranges z £ (0, 1), y £ (0, 00), x £ (0, 00). 
The index r > 1 takes in typical applications a half- 
integer value due to relativistic kinematics. It can be 
easily seen that d\f is a valid measure, positive definite, 
with bound integrals 



fi n — J dxW r (x; y, z)x n < 00 

for n a positive integer. As a consequence we can de- 
fine a family of orthogonal polynomials analogous to the 
Sonine polynomials, but more appropriate for a relativis- 
tic Bose-Einstein gas, which can conveniently be chosen 
monic (coefficient of highest dimension term equals one) , 
denoted P s (x; , y, z) and with an orthogonalization 



dxW r (x;y,z)ps(x;,y,z)P*(x;,y,z) = S st A s r (y,z) . 

(B2) 

Since the polynomials are considered monic, A is not 
unity. In the calculation presented we have considered 
a meson gas slightly out of equilibrium, where a good 
approximation is achieved by keeping only the first poly- 
nomial in the expansion of g{p) defined in i|30|) above. 
This has been verified in Q in the non-relativistic limit. 
Therefore in this calculation we only need to evaluate the 
Aoo(y,z) function. 



